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Abstract. Parallax measurements of pulsars allow for accurate measurements of the interstellar electron density 
and contribute to accurate tests of general relativity using binary systems. The Square Kilometre Array (SKA) 
will be an ideal instrument for measuring the parallax of pulsars, because it has a very high sensitivity, as well 
as baselines extending up to several thousands of kilometres. We performed simulations to estimate the number 
of pulsars for which the parallax can be measured with the SKA and the distance to which a parallax can be 
measured. We compare two different methods. The first method measures the parallax directly by utilising the 
long baselines of the SKA to form high angular resolution images. The second method uses the arrival times of the 
radio signals of pulsars to fit a transformation between time coordinates in the terrestrial frame and the comoving 
pulsar frame directly yielding the parallax. We find that with the first method a parallax with an accuracy of 
20% or less can be measured up to a maximum distance of 13 kpc, which would include 9 000 pulsars. By timing 
pulsars with the most stable arrival times for the radio emission, parallaxes can be measured for about 3 600 
millisecond pulsars up to a distance of 9 kpc with an accuracy of 20%. 
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1. Introduction 

Distance measurements have always played a hugely im- 
portant role in most aspects of astronomy, on both 
Galactic and extragalactic scales, but are often very dif- 
ficult to access observationally. Arguably, the most reli- 
able measurements are naturally obtained from parallax 
measurements when the annual movement of the Earth 
around the Sun can be used to detect a variation in the 
apparent source position. For nearly all sources and ap- 
plications, this involves imaging the source against the 
background sky. This can be done at optical wavelengths 
or, usually more precisely, at radio frequencies using Very 
Long Baseline Interferometry (VLBI). In the case of radio 
pulsars however, a second type of parallax can be mea- 
sured using timing measurements of the radio pulses. Here, 
distances are retrieved by detecting a variation in pulse ar- 
rival times at different positions of the Earth's orbit caused 
by the curvature of the incoming wavefront. In contrast to 
an imaging parallax, the precision of timing parallax mea- 



surements is highest for low ecliptic latitudes and lowest 
for the ecliptic pole. Pulsar timing of binary pulsars offers 
further possibilities to determine the distance from a sec- 
ular or annual variation in some orbital parameters. In all 
cases, however, the methods are usually limited to rela- 
tively nearby sources, but improvement in telescope base- 
line lengths and sensitivity promise improvements in pre- 
cision and hence the accurate measurement of distances. 

The most significant advance will be achieved by the 
Square Kilometre Array (SKA), which is being planned 
as a multi-purpose radio telescope constructed from many 
elements, both dishes and aperture arrays, leading to a 
total coll ecting area approach ing 1 million square me- 



tres (e.g. Dewdnev et al. 20091 ). In the frequency range 
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of 200 MHz to 3 GHz, the sensitivity of the SKA will be 
around 10000m 2 /K, but depending on the design that 
will be ch osen the sensitivity could be slightly higher or 
lower (see lSchilizzi et alJl2007h . The core of the SKA will 
be densely filled with 50% of the collecting area within a 
5 km diameter area. The remaining elements will be placed 
at locations extending up to several thousands of kilome- 
tres from the core. The field of view (FoV) of the dishes 
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will be around 0.64 deg 2 at 1.4 GHz for a single-pixel re- 
ceiver and might be extended up to 20 deg 2 by means of 
phased array feeds, although for the long baselines only 
the single-pixel receivers will be available. 

Having both a dense core and long baselines allows 
the SKA to both find radio pulsars and to perform ac- 
curate astrometric measurements on them. In particular, 
the SKA will be able to perform parallax measurements 
on a great number of pulsars, allowing accurate determi- 
nation of their distances from the Earth. Most distances to 
pulsars are currently derived from the measured column 
density of free electrons between the pulsars and Earth, 
known as the dispersion measure (DM). Using a model 
for the interstellar medium and its free electron content, 
the DM can be converted into a distance. However, the 
Galactic electron distribution is not accurately known, so 
that the distance derived from DM values is often uncer- 
tain, the error s in the distance es timates being as large as 

l2009h . 



100% 



see e.g. 



Delleretal 



The SKA will 



( Cordes et al. 200 



crform a Galactic census of pulsars 
that will potentially yield a popula- 
tion of 20,000 or more radio pulsars across the Galaxy 
(jSmits et al.l [20091. Parallax measurements of many of 



these pulsars with the SKA will provide direct measure- 
ments of their distances, which will then allow accu- 
rate studies of the ionised component of the interstellar 
medium, given the measurement of the DM. In this way, 
the interstellar electron model can be calibrated and will 
not only provide reliable estimates of the DM distances of 
pulsars without a parallax measurement in return, but will 
simultaneously provide a map of the free electron content 
in the Milky Way that can be combined with HI and HII 
measurements to unravel the Galactic structure and the 
distribution of ionised material. Combined with Faraday 
rotation measurements of the same pulsars, the Galactic 
magnetic field can also be studied in much greater detail 
as today, since the precise distance measurements poten- 
tially allow us to pin-point field-reversals occurring with 
some accuracy. 

The determination of accurate pulsar distances is also 
important for those pulsars that are part of binary sys- 
tems, in particular for thos e with another compact object 



such as the Double pulsar ([Burgav et al.ll2003t iLvne et al 
2004h . Relativistic effects can be used to determine the 



distance to some of these systems when the validity of 
general relativity is assumed. In reverse, to perform pre- 
cision tests of general relativity, kinematic effects have to 
be removed for which it is often required to know the 
distance precisely (lLorimer fc Kramerl2005 : Kramer et al 



20061 iDeller et al.ll2009h . 

In this paper, we review the possible ways to determine 
pulsar distances with the SKA. To compare the technical 
capabilities with the expected pulsar population, we sim - 



ulate an SKA pulsar survey, following lsmits et al.1 (l2009h . 



Using the results of this simulation we calculate the ac- 
curacy with which the SKA can determine the distance 
to each pulsar by measuring its parallax. We examine two 
basic methods for performing the astrometry, the first that 



uses imaging and the long baselines of the SKA, and the 
second that uses the timing of the pulsar radio signal. Here 
we distinguish between parallax measurements (applica- 
ble to all pulsars) and distance-related effects for binary 
pulsar parameters. 

In Sect. [U we describe the simulation of the SKA pul- 
sar survey. The two different methods for measuring the 
parallax of pulsars are explained in refsec:parallax, where 
we also explain how we estimate the capability of the SKA 
to measure the parallaxes of pulsars. The results are pre- 
sented in 21 Finally, we discuss our findings in[5j 

2. Simulation of an SKA pulsar population 

To demonstrate the SKA capabilities not only in princi- 
ple but to also estimate how many pulsars would offer the 
highest astrometry precision, we gauge the overall SKA 
performance for pulsar astrometry by simulating the pop- 
ulation of pulsars that can be expected to be involved in 
the experiments. 

The computational requirements for data analysis 
limit the amount of collecting area that can be used 
in an SKA p ulsar survey design ed to find new pulsars. 



According to ISmits et al.1 ((2009), at the beginning only 
the inner 1-km of the SKA core can be used, which con- 
tains ^20% of the total collectin g area. We per f ormed the 



same simulation as performed bv lLorimer et al. (2006), us- 
ing their Monte Car lo simulation package. In their study, 



Lorimer et al.l (|2006l ) used the results from recent surveys 



with the Parkes multi-beam system to derive an underly- 
ing population of 30 000 normal pulsars with an optimal 
set of probability density functions (PDFs) for pulsar pe- 
riod (P), 1400-MHz radio luminosity (L), Galactocentric 
radius (R), and height above the Galactic plane (z). 
For our simul a tions, we adopt the PDFs of model C in 



lLorimer et all (120061) . The intrinsic pulse width of each 



pulsar follows the power law relationsh ip with spin period 



given in Eq. (5) of lLorimer et al.l (|2006T) (see also their Fig 



4). To compute the expected DM and scatter broadening 
effects on each pulse, we use the NE2001 electron density 
model. To scale the scatter-broadening time to arbitrary 
SKA survey frequencies, we adop t a frequency power-law 
index of -4 (|Lohmer et al.ll200lh . 

We then simulated an all-sky pulsar survey at 800 MHz 
with a bandwidth of 500 MHz and 20% of the sensitivity 
of the full SKA (i.e. 2 000m 2 /K). We consider a pulsar 
to be detected by a given configuration if its flux density 
exceeds the threshold value, and its observed pulse width 
is less than the spin period. 

Since the core of the SKA will be located near the lat- 
itude of -30°, we limited the simulation to only detecting 
pulsars with a declination of 50° or less. This led to the 
detection of 15 000 normal pulsars in our simulation. In 
addition, the SKA would find 6 000 MSPs. Figure Q] shows 
the distribution of these pulsars in the Galaxy. The void 
at the centre of the figure is due to the extreme scattering 
of the radio signal at the Galactic centre, which decreases 
the sensitivity of detection. The degree of scattering near 
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Sgr A* is not yet known but we can expect that frequen- 
cies above 10 GHz may be required to detec t pulsars very 



close to the central supermassive black hole (jCordes et al 
20041) . 
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Fig. 1. Distribution of the 15 000 normal pulsars detected 
in the simulation (see text for details). The centre of the 
plot corresponds to the centre of the Galaxy. The grey 
areas mark the inner and outer disks and spiral arms used 
in the NE2001 electron density model. The black circle 
indicates the position of the Earth. 



3. Distance measurements 

The parallax and hence distance of a stellar object is clas- 
sically obtained by measuring the celestial coordinates of 
the object on the sky at different orbital positions of the 
Earth. Here we examine two methods for measuring the 
parallax of pulsars with the SKA, imaging parallax and 
timing parallax, before we also briefly discuss the "orbital" 
parallax measurements of binary pulsars. 

3.1. Imaging observations 

The straightforward method for measuring the parallax 
of pulsars is to measure the position of the pulsar on 
the sk y over time by means of imaging. Fomalont fc Reid 
(l2004h estimate the astrometric capabilities of the SKA by 
means of interferometry using the stated des ign goals of 
this proposed instrument ( Schilizzi et al.ll2~007t ). They out- 
line some of the calibration techniques for astrometry that 
could be employed by the SKA and discuss the effects of 
limited dynamic range on astrometry; they conclude that 
with multiple calibrator sources available at separations of 
a few arcminutes the astrometric accuracy that the SKA 
can potentially obtain is 9 a ~ 0/1000, where 9 a denotes 
the astrometric accuracy, 8 is the imaging resolution ob- 
tained by the array (typically FWHM of the synthesised 
beam) and the factor 1 000 is the dynamic range. At the 
frequency of 1.4 GHz and a 3 000 km maximum SKA base- 
line, the potential astrometric accuracy is approximately 



15/ias. However, for many pulsars the limiting factor will 
be given by the limited SNR of the pulsar detection. We 
therefore estimate the astrometric accuracy as the maxi- 
mum of 0/1000 and 6>/SNR pu i sar , where SNR pu i sar is the 
SNR of the pulsed flux from the pulsar over an integra- 
tion time of 12 hours. To achieve this SNR requires the 
correlator to perform de-dispersion and gating over the 
on-pulse of the pulsar. Therefore, an ephemeris of the pul- 
sar is required that can be obtained from regular timing 
observations. 

As of the end of 2009, parallaxes for approximately 30 
pulsars have been ob tained using VLBI (for a summary see 
IVerbiest et alj|2010l ). Typically about six measurements 
over a period of two years are required to measure the 
parallax and proper motion of a single pulsar, typically 
requiring a total observing time of 72 hours of observa- 
tional time per pulsar (12 hours per VLBI observation). 
Measuring the parallax of multiple pulsars in the FoV of 
the telescope involves using the correlator to perform de- 
dispersion and gating for every pulsar in the FoV, which 
is possible for current VLBI observations using the capa- 
bilities of new software correlation systems. 

3.2. Timing observations 

A second set of methods for performing astrometry of a 
pulsar involves the accurate timing of the pulsar radio 
signals, i.e. the measurement of the pulse times-of-arrivals 
(TOAs) at the telescope. Here we can distinguish between 
three methods, all relying on very precise TOA determina- 
tion: a) parallax measurements using the Earth orbit, b) 
parallax measurements for binary pulsars using the Earth 
orbit and that of the pulsar, and c) distance estimates of 
binary pulsars based on the comparison of observed orbital 
parameters with those predicted by general relativity. 



3.2.1. Timing parallax 

The "classical" timing parallax measurement utilises the 
fact that the wave front curvature of a pulsar signal is 
directly related to the distance of the source. For an in- 
finitely distant source, the wave front is planar and, after 
taking into account Earth's orbit, the TOA would essen- 
tially be the same at all orbital phases of the Earth. For 
a finite distance, however, the curvature of the wavefront 
introduces an annual periodic change in the apparent di- 
rection, hence a six-monthly periodicity in the TOAs due 
to the changing path - length from the pulsar to the Earth 
(|Lvne fc Smith! 12004 iLorimer fc Kramerl [2005h . The ap- 
parent change in direction is more easily measurable for 
low ecliptic latitudes - in contrast to an imaging paral- 
lax - as the amplitude of the TOA variation scales with 
cos 2 (3/d, where /3 is the ecliptic latitude and d the distance 
to the pulsar. 

Performing simulations of parallax measurements (see 
Fig 12]) deviations from the simple cos 2 (3 scaling become 
apparent: it is still possible to measure a parallax with 
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Fig. 2. Left) Post-fit residuals for a pulsar at a distance of 200 pc located at the pole of the ecliptic after fitting for 
position and proper motion but not for parallax. Remaining structure in the TOAs is clearly visible. Right) Error in 
the timing parallax as a function of the ecliptic latitude /? for weekly timing of a pulsar to a TOA with 10 ns error. 
The blue and red lines correspond to 5 years and 10 years of timing, respectively. 



a finite precision even at high ecliptic latitudes. As the 
Earth's orbit deviates from a circular shape because of a 
small eccentricity and deviations caused by other masses, 
non-linear terms depending on e 2 cause differences in the 
arrival times that are measurable if the precision is suffi- 
ciently high. 

Figure [5] shows the absolute uncertainty in a paral- 
lax measurement for an observing span of 5 years and 10 
years, respectively. In both cases, weekly TOAs with 10 ns 
precision are assumed, which is challenging but not impos- 
sible with SKA sensitivity (see Section [5]). In these cases, 
it will be possible to determine the parallax of a pulsar 
with a precision of ~ 1/^as for pulsars near the ecliptic, 
and still with a precision of 300 fi&s and 70 /ias, respec- 
tively, at the ecliptic pole. In other words, in principle the 
SKA should allow us to achieve a distance measurement 
with an uncertainty of less than 5% out to about 30 kpc! 
If the TOA errors are larger, we can still obtain a pre- 
cision of 20% out to 30 kpc with 50 ns TOA errors and 
to 15 kpc with 100 ns TOA errors, respectively. While a 
timing precision of 100 ns has been achieved already for 
some pulsars, in reality the timing precision will decrease 
with distance as the pulsars become fainter. We discuss 
this further in Sect. |U 

3.2.2. Orbital parallax 

If a pulsar happens to be in a binary system, the pulsar 
orbit will be viewed under slightly different angles from 
different positions of the Earth's annual orbit. The result 
is a periodic change in the observed longitude of peri- 
astron, w, and the projected semi-major axis x = asinz, 
where a is the semi-major axis of the puls ar orbit. This ef- 
fect, known as the annual o rbital parallax (IKopeikinlll995 ; 
van Straten fc Bailesl l2003f) , causes a variation in x and to 
given by 



w(i) 



UJQ 



-rssB 



I' 



(2) 



x(t) 



x-a 



coti 



TSSB (t) " J' 



(1) 



where J' = — sin57 asc io + cos£l asc .Jo and I' = 
— cosf2 asc Jo + sinf2 asc Jo- The vectors Jo and Jq and the 
longitude of the ascending node f2 asc , d efine the orbital ori- 
entat ion in space (see e.g. Figure 8.3 in lLorimer fc Kramer 
120051 ). Since the amplitude of this effect depends on the 
distance d, a detection of the annual orbital parallax can 
also lead to an accurate pulsar distance measurement if 
the orientation of the orbit can be determined, e.g. from 
aberration and/or spin-precession effects. Furthermore, a 
significant proper motion in the sky can lead to measur- 
able secular changes in x and w, which ca n be used to con - 
strain the orientation of the pulsar orbit (|Kopeikinlll996l ) . 
For this, one needs to ensure that the other contributions 
to these secular changes are either well understood or well 
below the measurement precision of x and w. For instance, 
the secular variation in x can be the result of various ef- 
fects that can be summarised as: 



PM 



, GW , 
X \ d£A 



D 
D' 



(3) 



where the contributions to the observed semi-major axis 
derivative are a change in the orbital inclination due to the 
proper motion, as outlined above, the emission of gravita- 
tional waves (GW), a varying aberration, d£A/di, and a 
changing (radial) velocity contribution, —D/D. 

The last term of Eq. ([3]) is a combined effect of a chang- 
ing Doppler shift caused by a relative acceleration in the 
gravitational field of the Galaxy or a globular cluster and 
an appa rent acceleration due to the proper motion across 
the sky (jDamour fc Tavlorlll99ll) 



D c 



(dPSR — a SSB) 



(4) 



where Kq is the unit vector from the Earth towards the 
pulsar, and apsR and ossb are the Galactic accelerations 
at the location of the binary system and the Solar System 
barycentre. The last term including the transverse veloc- 
ity, Vt, is known as the Shklovskii term. Consequently 
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D/D depends on the distance to the system, and in prin- 
ciple could provide another way to estimate the distance 
to the pulsar. In practice, this term is generally below the 
measurment precision, and anyway difficult to separate 
from the other contributions. However, the acceleration 
D/D can play an important role in the distance measur- 
ment as outlined in the next section. 



le-13 
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Fig. 3. Decrease in orbital period for a pulsar in a 1-d 
orbit with a 0.3Mq companion and a transverse velocity of 
100 km s _1 . The green curve shows the contribution from 
gravitational wave damping, while the blue curve is drawn 
for different time spans of observations, namely 5, 10, and 
15 yr, respectively. We assume weekly observations, and a 
10 ns TOA uncertainty that increases quadratically with 
distance to account for the decreasing flux density. 



where all masses are expressed in solar units by using 
the constant T© = (GM Q /c 3 ) = 4.925490947 fjs, G is 
Newton's gravitational constant, c the speed of light, and 



/(e) 



[1 + (73/24)e 2 + (37/96)e 4 
(1 -e 2 ) 7 /2 



(6) 



However, the observed value of Pb will often differ from 
the expected intrinsic value (Eq. [S]) by an additional term 
that is again given by the relative acceleration, namely 




obs 




GW 



(7) 



where we have ignored possible contri butions due to mass 
loss o r tidal effects in the system (see lLorimer fc Kramer 
20051) . 
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3.2.3. Using general relativity for distance 
measurements or kinematic parallax 

While the observed x can be contaminated by a variety of 
effects (see Eq. [3]), the observed orbital decay rate of an 
orbit, Pb, is not affected by aberration or proper motion 
effects. If the intrinsic decay rate is determined purely by 
the loss of orbital energy due to GW emission, one can 
predict Pb if the pulsar and companion masses are ob- 
tained via pulsar timing thanks to the measurement of 
post-Keplerian (PK) parameters. If two PK parameters 
are inferred, for instance from relativistic orbital preces- 
sion, gravitational redshift, or a Shapiro delay, the m asses 
can be determined (see e.g. lLorimer fc Kramerll2005l )) and 
we can compute the value of Pb expected for GW emission 
assu ming that general relativity (GR ) is correct. In this 
case ( Damour fc Deruelldll985lll986t l. 



pGW 



1927T 5/3 

i- c\ 



2^ /(6) 



(m F 
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Fig. 4. Precision for distance measurements inferred from 
the timing of a pulsar at 1 kpc distance as a function of 
observing time span. The red curve is for the "classical" 
timing parallax at two different ecliptic latitudes, j3. The 
distance uncertainties achieved with the kinematic paral- 
lax are shown in blue for two different orbital periods. We 
again assume weekly observations with a TOA precision 
of 10 ns. 

The term D/D can be determined from Eq. [7] by com- 
paring the computed GR value with the observations. 
Figure [3] compares the relative size of the contributing 
terms by showing the orbital period change of a pulsar in 
a 1-d orbit with a O.3M companion and a typical trans- 
verse velocity of 100 km s -1 . The red curve shows the 
kinetic contribution of the Shklovskii term, the green one 
representing the gravitational wave damping term and the 
blue line compares the timing precision of three different 
time spans of observations, namely 5, 10, and 15 years, 
respectively. Figure [4] compares the resulting absolute un- 
certainty in distance measurement for the "classical" tim- 
ing parallax (red) at two different ecliptic latitudes with 
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the results for the "kinematic" parallax (blue) with dif- 
ferent orbital periods. It is interesting to see that for the 
first few years of observations the classical timing parallax 
yields much higher precision but that in the long-term the 
kinematic parallax eventually leads to more precise results 
assuming that the acceleration in the Galactic gravita- 
tional potential can be accounted for and that the masses 
are known with sufficient accuracy via the measurement 
of PK parameters. 

4. Results 

While we have demonstrated the SKA capabilities in prin- 
ciple, we now review the applicability of our results to real 
SKA observations. 
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plane and 450 days to measure the parallaxes of all the 
detectable pulsars to a parallax error of 20%. 




ji/Aji 



Fig. 5. Histogram of ir/ Air for the pulsars detected in the 
simulation of imaging observations. The vertical dotted 
line marks the tt/Att — 5 cutoff. 



4.1. Imaging observations 

Fig.[5]shows the histogram of the quantity tt/Att, where tt 
is the parallax and Att is the estimated error in the paral- 
lax, for the simulated sample. The histogram corresponds 
to two underlying distributions for which the parallax is 
limited by the SNR of the pulsar or the dynamic range. 
The low number of pulsars around tt/Att ~ 1 arises from 
the detection limit of the SKA survey. The main peak cor- 
responds to the large number of pulsars with a low SNR. 
The rise in the distribution at tt/Att ~ 5 reflects the peak 
of the pulsar distance distribution at ~6 kpc. Taking a 
cutoff tt/Att = 5, we find that the SKA can potentially 
measure the parallaxes for ~ 9,000 pulsars with an error 
of 20% or smaller. This includes pulsars out to a distance 
of 13 kpc. For the Galactic plane, as many as 25 pulsars 
can be observed in the FoV when using a single-pixel re- 
ceiver. If we assume that the correlator can deal with all 
the pulsars in the FoV, it will take 215 days to measure 
the parallaxes of all the pulsars (9 000) in the Galactic 
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Fig. 6. Fractional precision of distance measurement by 
timing parallax for weekly observations. The TOA error 
is assumed to be 10 ns up to 3 kpc beyond which the error 
scales with distance squared. The blue line indicates five 
years of observations, the red one ten years, respectively. 



4.2. Timing observations 

We have shown that timing observations of pulsars can 
yield a parallax precision that exceeds that of an imag- 
ing parallax by at least an order of magnitude in terms 
of both relative and absolute precision. For longer time 
spans, the kinematic parallax is of even higher accuracy. 
For these results, we have typically assumed a TOA un- 
certainty of 10 ns. Such a precision is only possible for 
recycled millisecond pulsars (MSPs) and is a factor of 
5-10 higher than what has been achieved for the best 
MSPs to date. Increased sensitivity is expected to help 
improve the timing precision in a linear way, provid- 
ing an up to 100 times higher sensitivity with the SKA 
and ensuri ng that such timing prec ision appears easily 
achievable ([Lorimer fc Kramerl 120051 ) . Uncertain factors 
are, however, possible intrinsic timing noise of the pulsar 
clock, propagation effects in variable interstellar weather, 
and both instrumental and improper calibration effects. 
Recent and ongoing studies ( Verbiest et al. 20091 Kuo et 
al. in prep.) indicate that despite challenges, a timing pre- 
cision for MSPs between 10 and 50 ns seems feasible. We 
expect that about ~ 6000 of the pulsa rs detected with the 
SKA will be MSPs (jSmits et alJl2009h but not all of those 
will be suitable for high precision timing. Moreover, as the 
new pulsars become more distant, the flux density will de- 
crease and hence the timing precision becomes worse. We 
attempt to take this into account by increasing the TOA 
errors by a factor of (d/do) 2 , where do determines the 
distance beyond which the timing precision is dominated 
by the radio meter equation (rather than other limiting 
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factors). Inspecting the timing precision of the presently 
known MSPs as a function of known distance, we estimate 
this transition to be around 3 kpc. This number may be 
significantly larger for the SKA, so that the present re- 
sults using classical timing parallax shown in Fig. [5] are 
probably a conservative estimate. These suggest that dis- 
tances can be measured with a precision of at least 10% 
to a distance of about 10 kpc. 
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Fig. 7. Fractional precision of a distance measurement for 
weekly observations with TO A errors of 10 ns up to a 
distance of 3 kpc. Beyond 3 kpc, we assume the timing 
precision to become worse by a factor scaling with the 
square of the distance. Colours indicate the length of the 
observing span: 5 yr (blue), 10 yr (red), and 15 yr (green), 
respectively. The thick lines indicate the kinematic dis- 
tance measurement for a pulsar in a 1 day orbit with an 
0.3 solar mass companion. The thin lines indicate the tim- 
ing parallax for the same pulsar. Note that the 10-yr lines 
(red) overlap. 

In Fig. [71 we compare the distance measurements using 
classical timing parallax and using kinematic parallax and 
adjusting the TOA precision for the distance scaling. As 
before, the classical timing parallax yields the superior 
results for the first years of observations, but is overtaken 
by the kinematic parallax after 10 years. 

To gauge how many MSPs would be likely to reach 
the timing precision assumed above, we compare the ob- 
tained precision among a consistent group of known MSPs 
at a relatively well-determined distance by inspecting the 
timing results for those 16 MSPs in the globular cluster 
47 Tucanae that are regularly observed. The other eight 
MSPs in the cluster are at the current limit of detectabil- 
ity of the Parkes telescope and can only be timed when 
the interstellar scintillat ion boosts the flu x density above 
the detection threshold (jFreire et alJl2003l ). Tabled] shows 
the current and predicted TOA uncertainty as well as an 
estimate of the accuracy with which the parallax of these 
16 MSPs can be measured with the SKA. These estimates 
are based on a scaling that ensures that the TOAs are 



improved by a factor of ~ 100. However, as the full SKA 
may not be able to be phased up for timing observations, 
we conservatively estimate that we can only use the inner 
5 km, containing 50% of the collecting area, so that we 
scale the TOA errors by a factor of 50. To estimate the 
parallax error inferred from timing each of the 6 000 new 
MSPs that will be found by the SKA, we randomly as- 
sign a TOA error from the 16 MSPs in 47 Tucanae. The 
parallax error is then corrected for ecliptic latitude and 
distance following the results shown in Figs. [5] and [7J 
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Fig. 8. Histogram of ir/Air for the pulsars detected in the 
simulation of timing 6 000 MSPs. The vertical dotted line 
indicates the ir/Air = 5 cutoff. 

Figure \E\ shows the histogram for the quantity tt/Att, 
where tt is the parallax determined by means of timing 
MSPs and An is the estimated error in the parallax for the 
simulated sample. Taking once again a cutoff tt/Att = 5, 
we find that the SKA can potentially measure the paral- 
laxes for ~ 3 600 MSPs with an error of 20% or smaller out 
to a distance of 9 kpc. Timing all these MSPs using the 
phased array feeds of the SKA with an observation time 
of 1 hour, will take roughly 25 days. To obtain a timing 
solution requires re-observations of these MSPs once ev- 
ery two weeks over a period of at least one year. The total 
observation time would then amount to 650 days. If we 
were to time these MSPs only to achieve a parallax error 
of 20%, the total observation time would drop to roughly 
40 days. In any case, there will be a follow-up timing for 
each MSP that the SKA will detect to select the most in- 
teresting MSPs. Therefore, obtaining parallaxes for many 
of these MSPs will not require extra observation time. 



4.3. Pulsar proper motions 

Astrometry with the SKA will permit us to obtain accu- 
rate measurements of the proper motion of pulsars. These 
are extremely valuable for constraining the origin of pul- 
sar velocities, which is currently known with an accuracy 
of only a few to tens of mas per year for normal pulsars 
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Table 1. Error in TOAs from 16 MSPs from the cluster 
47 Tucanae and the resulting error in parallax from the 
timing solution. 
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6 


0.12 
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9 


0.18 


31 


J0024- 
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10 


0.20 


34 


J0024- 


7204G 


7 


0.14 


24 


J0024- 


7204H 


20 


0.40 


69 


J0024- 


72041 


14 


0.28 


88 


J0024- 


7204J 


3 


0.06 


10 


J0024- 


7204L 


32 


0.64 


110 


J0024- 


7204M 


17 


0.34 


58 


J0024- 


7204N 


12 


0.24 


41 


J0024- 


7204O 


8 


0.16 


28 


J0024- 


7204Q 


19 


0.38 


65 


J0024- 


7204S 


11 


0.22 


38 


J0024- 


7204T 


49 


0.98 


168 


J0024- 


7204U 


9 


0.18 


31 



and o ften < 1 mas per year for MSPs (jLorimer fc Kramer 
2005). Using imaging observations with the SKA allows 
a potential accuracy of 15 fj,as per year. From our simula- 
tions, we find that for 60% of all pulsars the proper motion 
can be measured with an accuracy of 15 — 25 /xas per year. 
Furthermore, we find that using timing observations, the 
SKA can achieve a proper motion accuracy of < 5/ias per 
year for almost 1 000 MSPs and an accuracy of < 30/ias 
per year for more than 3 000 MSPs. 



5. Discussion 

We have investigated two methods to measure the paral- 
lax of radio pulsars using the SKA. The imaging parallax 
method utilises the long baselines of the SKA to keep track 
of the position of the pulsar on the sky during different 
orbital positions of the Earth. The timing method uses 
the arrival times of the radio signals of pulsars to derive 
a "classical" parallax, an "orbital-annual" parallax, or a 
"kinematic" parallax. With imaging parallax, it is possible 
to determine parallaxes with an error of 20% or smaller 
for 9 000 pulsars out to a distance of 13kpc. Measuring 
the parallaxes of all the pulsars in the Galactic plane 
alone will require 215 full days of observing. Furthermore, 
to obtain the required astrometric accuracy of 15^as re- 
quires several calibrators to be within a few arc-minutes of 
the source. With the timing parallax method, parallaxes 
can be measured out to 9 kpc for about 3 600 MSPs, with 
20% or higher accuracy. This would require a minimum of 
roughly 40 days of observation time. However, there will 
have to be follow-up timing in any case for each MSP that 
will be detected by the SKA to select the most interesting 



MSPs. Therefore, obtaining parallaxes for many of these 
MSPs will not require extra observation time. 

When we consider the increased timing precision with 
the SKA and the recent realization that the apparent 
"timing noise" in non-millisecond pulsars is non-random 
and o f magnetospheric origin, possibly allowing to correct 
for it (jLvne et al.ll2010l ) , it seems likely that parallax mea- 
surements will even be possible for normal, non-recycled 
pulsars. In this case, the sample of available pulsars could 
easily exceed the number of 10 000 but this is difficult to 
estimate at this stage. 

The imaging and timing parallax methods are com- 
plementary, with timing parallax providing high quality 
measurements for a relatively small number of objects but 
to large distances, whereas imaging parallax can provide 
good results for a large number of closer objects. Timing 
parallax is limited to millisecond pulsars with stable tim- 
ing characteristics, whereas an imaging parallax depends 
only on the observed strength of the radio emission, not 
the rotation characteristics of the pulsar. Imaging paral- 
lax measurements can indeed help us significantly with 
precision tests of general relativity to correct for the same 
kinematic effects that we can otherwise use for precise 
distance measurements. Moreover, obtaining imaging in- 
formation immediately after the discovery of a new pulsar 
will help us to obtain positional information that will help 
the typical "solving procedure" of a pulsar, i.e. the first 
determination of a coherent timing solution. In general, 
solving pulsars takes several weeks to months of intensive 
observing, but the positional information from imaging 
eradicates the correlations between position and period 
derivative and allows us to reduce the cadence of the ob- 
servations significantly. 

From the above results, it is clear that the SKA will be- 
come a superb astrometry instrument, the data of which 
feeds directly back into fundamental astrophysical ques- 
tions. For instance, the distance to about 9 000 pulsars 
can be measured with an error of 20% or smaller. The 
dispersion measure from each pulsar can then be used 
to accurately calculate the electron density between the 
Earth and each of these pulsars. This will result in a map 
of the distribution of the ionised gas of high accuracy out 
to a distance of 13 kpc. By timing the most stable MSPs 
that the SKA will discover, the ionised gas can be mapped 
to a distance of 9 kpc. By imaging pulsars, we can locate 
them with a precision that is comparable with that of 
GAIA, which is able to measure distances out to 10 kpc 
with a precision of 10%. This is also achievable by tim- 
ing MSPs and may be able to be, for selected sources, 
improved upon. This is exciting as we will be able to get 
independent distance measurements from GAIA and the 
SKA for pulsars that have an optical companion or are 
located, e.g. in a globular cluster. 

The simulations presented in this paper are no t cor- 
rected for the Lutz-Kelker bias (jLutz fc Kelkerl 19731) . This 
bias is introduced when the parallax measurement does 
not take into account the larger volume of space that is 
sampled at smaller parallax values and leads to a sys- 
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tematic overestimate of the parallax, thus to an underes- 
timate of the distance to the star. This bias is strongly 
dependent on the accuracy of the parallax measurement, 
becomin g smaller as th e meas urement becomes more ac- 
curate. IVerbiest et ail (|2010h have studied this bias for 
the specific case of the parallax measurements of neutron 
stars, incorporating the bias introduced by the intrinsic 
radio luminosity function and a realistic Galactic popu- 
lation model for neutron stars. They found a significant 
bias for measurements with a 50% error. In the present 
paper, we consider parallax accuracies of 20% and higher. 
We find that in our simulations the Lutz-Kclker bias is 
always smaller than one standard deviation and can be 
neglected for most pulsars, hence does not affect our con- 
clusions. Nevertheless, since the bias in a pulsar parallax 
measurement depends on the luminosity of the pulsar and 
the position on the sky, it is advisable to correct for it in 
every pulsar parallax measurement^. 

Overall, our result shows that the SKA requires both 
a concentration of collecting area in the core, to search 
for and detect the Galactic pulsar population and then 
measure accurate timing parallaxes, and a distribution of 
collecting area on long baselines, to achieve the angular 
resolution required to perform imaging parallax measure- 
ments. Thus, the current design specifications for the SKA 
are adequate for achieving both of these objectives. 



Kopeikin, S. M. 1996, ApJ, 467, L93+ 

Kramer, M., Stairs, I. H., Manchester, R. N., et al. 2006, sub- 
mitted 

Lohmer, O., Kramer, M., Mitra, D., Lorimer, D. R., & Lyne, 

A. G. 2001, 562, L157 
Lorimer, D. R., Faulkner, A. J., Lyne, A. G., et al. 2006, 

MNRAS, 372, 777 
Lorimer, D. R. & Kramer, M. 2005, Handbook of Pulsar 

Astronomy (Cambridge University Press) 
Lutz, T. E. & Kelker, D. H. 1973, PASP, 85, 573 
Lyne, A., Hobbs, G., Kramer, M., Stairs, I., & Stappers, B. 

2010, Science, 329, 408 
Lyne, A. G., Burgay, M., Kramer, M., et al. 2004, Science, 303, 

1153 

Lyne, A. G. & Smith, F. G. 2004, Pulsar Astronomy, 3rd ed. 

(Cambridge: Cambridge University Press) 
Schilizzi, R. T., Alexander, P., Cordes, 

J. M., et al. 2007, SKA draft, 

http : //www. skat elescope . org/PDF/Preliminary_SKA_Specif icatio: 
Smits, R., Kramer, M., Stappers, B., et al. 2009, A&A, 493, 

1161 

van Straten, W. & Bailes, M. 2003, in Astronomical Society of 

the Pacific Conference Series, Vol. 302, Radio Pulsars, ed. 

M. Bailes, D. J. Nice, & S. E. Thorsett, 65-68 
Verbiest, J. P. W., Bailes, M., Coles, W. A., et al. 2009, 

MNRAS, 400, 951 
Verbiest, J. P. W., Lorimer, D. R., & McLaughlin, M. A. 2010, 

MNRAS, 405, 564 



Acknowledgements. The authors would like to thank an anony- 
mous referee for his/her comments. This effort /activity is sup- 
ported by the European Community Framework Programme 
6, Square Kilometre Array Design Studies (SKADS), con- 
tract no 011938. We gratefully acknowledge support from 
ERC Advanced Grant "LEAP", Grant Agreement Number 
227947 (PI Michael Kramer). The International Centre for 
Radio Astronomy Research is a Joint Venture between Curtin 
University and The University of Western Australia, funded by 
the Western Australian State Government. 



References 

Burgay, M., D'Amico, N., Possenti, A., et al. 2003, Nature, 
426, 531 

Cordes, J. M., Kramer, M., Lazio, T. J. W., et al. 2004, 48, 
1413 

Damour, T. & Deruelle, N. 1985, Ann. Inst. H. Poincare 

(Physique Theorique), 43, 107 
Damour, T. & Deruelle, N. 1986, Ann. Inst. H. Poincare 

(Physique Theorique), 44, 263 
Damour, T. & Taylor, J. H. 1991, ApJ, 366, 501 
Deller, A. T., Bailes, M., & Tingay, S. J. 2009, ArXiv e-prints 
Dewdney, P. E., Hall, P. J., Schilizzi, R. T., & Lazio, T. J. L. W. 

2009, IEEE Proceedings, 97, 1482 
Fomalont, E. & Reid, M. 2004, New Astronomy Review, 48, 

1473 

Freire, P. C, Camilo, F., Kramer, M., et al. 2003, MNRAS, 

340, 1359 
Kopeikin, S. M. 1995, 439, L5 



1 The Lutz-Kelker bias for pulsars can be determined online: 
http : //psrpop . phys . wvu . edu/LKbias 



